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ABSTRACT 



O 

O ; The NRAO VLA and VLBA were used from 1988 November to 1998 June to 

CNJ ' monitor the radio continuum counterpart to the optical broad line region (BLR) in the 

^^ ■ Seyfert galaxy NGC 5548. Photometric and astrometric observations were obtained 

at 21 epochs. The radio nucleus appeared resolved, so comparisons were limited to 
J^ . observations spanning 10-60 days and 3-4 yr, and acquired at matched resolutions of 

CNl . 1210 mas = 640 pc (9 VLA observations), 500 mas = 260 pc (9 VLA observations), or 

^r I 2.3 mas = 1.2 pc (3 VLBA observations). The nucleus is photometrically variable at 

8.4 GHz by 33 it 5% and 52 it 5% between VLA observations separated by 41 days and 
S»^ ' 4.1 yr, respectively. The 41-day changes are milder (19 it 5%) at 4.9 GHz and exhibit 

■^ ' an inverted spectrum (a ~ +0.3 it 0.1, S oc v^") from 4.9 to 8.4 GHz. The nucleus 

Q^' is astrometrically stable at 8.4 GHz, to an accuracy of 28 mas = 15 pc between VLA 

Q . observations separated by 4.1 yr and to an accuracy of 1.8 mas = 0.95 pc between 

J^ , VLBA observations separated by 3.1 yr. Such photometric variability and astrometric 

j^ I stability is consistent with a black hole being the ultimate energy source for the BLR, 

but is problematic for star cluster models that treat the BLR as a compact supernova 
k^ ^ remnant and, for NGC 5548, require a new supernova event every 1.7 yr within an 

Vh ' effective radius r^. = 42 mas = 22 pc. A deep image at 8.4 GHz with resolution 



660 mas = 350 pc was obtained by adding data from quiescent VLA observations. This 
image shows faint bipolar lobes straddling the radio nucleus and spanning 12 arcsec = 
6.4 kpc. These synchrotron-emitting lobes could be driven by twin jets or a bipolar 
wind from the Seyfert 1 nucleus. 

Subject headings: galaxies: individual (NGC 5548, UGC 09249) - galaxies: nuclei - 
galaxies: Seyfert - radio continuum: galaxies 



1. Motivation 

The Seyfert 1 galaxy NGC 5548 was the target of intensive space- and ground-based 
monitoring designed to measure the light travel time across its broad optical-line-emitting region 



(Peterson 1993, Korista et al. 19"95| , and references therein). These reverberation studies imply 
that the broad hne region (BLR) has a characteristic diameter ^ 40 light days, or < 64 //as 
for Hq = 50 km s~^ Mpc~^ and the galaxy's recession velocity relative to the 3 K background 
(5354 km s~^; de Vaucouleurs et al. 1991] ) . Given this small size, is it possible to constrain models 
for the ultimate energy source for the BLR in NGC 5548? The common view is that the energy 
source is a supermassive black hole fuelled by an accretion disk (e.g., Peterson 1997] ), and three 
recent findings support this scenario for NGC 5548. First, analysis of variations in the emission 
line profiles of NGC 5548 suggests that the BLR material moves along randomly inclined Keplerian 
orbits defined by a central mass of order ~ 7 x 10^ Mq ( [Wanders et al. 1995 , Wanders Sz Peterson 



1996, Peterson & Wandel 1999). Second, the central continuum source must irradiate the BLR 



anisotroptically ( Wanders et al. 1995| , [Wanders Sz Peterson 199^ ), perhaps because the continuum 
source is surrounded by an optically-thick accretion disk (cf. Figure 1 of poad &: Wanders 1996 ). 
Third, NGC 5548 exhibits very broad and asymmetric X-ray emission lines whose profiles are well 



matched by gravitationally-redshifted and fluorescing lines from accretion disk material (Nandra 
et al. 1997[ ). 

Still, the starburst hypothesis of Terlevich et al. (1992) potentially offers an alternative energy 
source: a compact star cluster hosts a supernova and the BLR is identifled with the compact 
supernova remnant formed as the supernova ejecta interacts with dense circumstellar material. 
Although the star cluster hypothesis has yet to be tested against the three findings listed above 
for NGC 5548, it appears consistent with analysis of the optical continuum variability and the 
reverberation results ( Aretxaga fc Terlevich 1993 , [Terlevich et al. 1995[) . Astrometric monitoring 
of Seyfert nuclei offers another strong test of of the starburst hypothesis ( [Melnick, Aretxaga, fc 
[Terlevich 1997 ). The low supernova rates and their short radiative lifetimes mean that as the 
compact star cluster evolves and new supernovae occur, the BLR should dance within the star 
cluster. In constrast, black hole models for low-power objects like Seyfert nuclei predict no motion 
of the nucleus over time, or simple linear motion if the emission arises from propogating shocks in 
a jet and counterjet ( Falcke fc Biermann 1996 , Yi &: Boughn 1998 , Yi fc Boughn 199S| ). This paper 
imposes observational constraints on the astrometric stability of the radio continuum nucleus of 
NGC 5548, to help distinguish between star cluster and black hole models for the energy source of 
the BLR. Although the star cluster models do not predict light curves for the radio continuum 
emission from a compact supernova remnant, Terlevich et al. (1995) do remark on the relevance 
of SN 1988Z to their models, as being an example of a (radio) supernova expanding into dense 
circumstellar material ( ptathakis fc Sadler 1991 , Van Dyk et al. 1993| ). Therefore this paper also 
examines the photometric variability of the radio continuum nucleus of NGC 5548, for comparison 
with the photometric evolution of SN 1988Z at radio wavelengths. 



2. Astrometric Imaging 

Table 1 gives a log of the astrometric imaging observations of NGC 5548 with the Very Large 



Array (VLA, Thompson et al. 1980) and the Very Long Baseline Array (VLBA, Napier et al. 
1994 ). Both instruments acquired data in dual circular polarizations at the center frequencies 
tabulated. For each polarization, the VLA bandwidth was 100 MHz with 3-level sampling and the 
VLBA bandwidth was 32 MHz with 4-level sampling. All VLA and VLBA observations after 1990 
assumed a coordinate equinox of 2000. Earlier VLA observations assumed an equinox of 1950 and 
were transformed to 2000 using a reference date of 1979.9 for the phase calibrator grid used at the 



VLA prior to 1990 ( [Walker 1999| ). 



On the VLA, an 8-minute observation of NGC 5548 was preceded and followed by an 
observation of the phase cahbrator OQ208 (B1404+286) about 4° from NGC 5548 (Figure 1). 
VLA observations of 3C286 were used to set the amplitude scale to an accuracy of about 3%. On 
the VLBA, a 3-minute observation of NGC 5548 was preceded and followed by an observation of 
the phase, rate, and delay calibrator J1419+27 about 2° from NGC 5548 (Figure 1). The position 
assumed for J1419+27 was aj2ooo = 14''19™59!29698 and 5j2ooo = +2 7° 06' 2 5'.' 5 540, with an error 
in each coordinate of less than 1 mas (M. Eubanks, private communication). Sources OQ208 and 
J1310+32 were also observed with the VLBA to check the astrometric accuracy and to align the 
phases of the baseband channels, respectively. VLBA system temperatures and gains were used 
to set the amplitude scale to an accuracy of about 5%, after first correcting for sampler errors. 
For both the VLA and VLBA, the position assumed for OQ 208 was aj2ooo = 14''07'^00!3943640 
and (5j2ooo = +28°27'14'.'69129 (M. Eubanks, private communication). These assumed positions 
for both J1419+27 (B1417+273) and OQ208 agree with those in Ma et al. (1998). The VLA data 
were calibrated using the 1994 January 15 release of the NRAO AIPS software, while the VLBA 
data were calibrated with the 1999 April 15 release. 

The AIPS task IMAGR was used to form and deconvolve images of the Stokes / emission 
from NGC 5548. Neither the VLA nor VLBA data on NGC 5548 were self-calibrated. The VLA 
images were made with a weighting scheme that achieved high angular resolution without seriously 
degrading sensitivity. Furthermore, VLA imaging tests showed that the peak source strength 
increased as angular resolution decreased, suggesting that the source was slightly resolved. Thus 
only VLA images with matched resolutions were analysed. The VLBA images were made with 
natural weighting to optimize sensitivity and were restored with a matched elliptical Gaussian 
beam. Photometric and astrometric comparisons will be limited to observations at matched 
angular resolutions of 500 mas = 260 pc (Figure 2, 9 VLA observations), 1210 mas = 640 pc 
(Figure 3, 9 VLA observations), and 2.3 mas = 1.2 pc (Figure 4, 3 VLBA observations). 

Image rms values (o"/) appear in Table 1 in units of microjanskys per beam area (^Jy ba~^); 
all rms values match those expected given the array and observation parameters ( |Wrobel fc Walker 



19991) . The image resolutions, expressed either as the FWHM dimension of a circular Gaussian 



beam or the FWHM dimensions and elongation orientation of an elliptical Gaussian beam, also 



appear in Table 1. The tabulated position peaks were derived from quadratic fits in the image 
plane. For the VLA images, these quadratic fits also yielded the peak flux densities listed in the 
table, along with errors that are quadratic sums of a 3% scale error and aj. For the VLBA images, 
the source appeared to be resolved, so total flux densities were derived by integration over A^ ~ 10 
beam areas and tabulated along with errors that are quadratic sums of a 5% scale error and 
aj\/N. This apparent resolution could be due either to unmodelled atmospheric and/or station 



motion effects ( Beasley &: Conway 1995 ), or due to real structure. That the structure is similar in 
the VLBA images suggests the latter but does not prove it, because some phase-referencing errors 
can be systematic from observation to observation. 

The table also notes the laA error in each coordinate for the astrometry of the peak position. 
a A is the quadratic sum of a term measuring the peak signal-to- noise ratio ( [Ball 1975 ) and a 



term measuring the accuracy of the absolute astrometry. This latter term was estimated in the 
VLA case as follows. The imaging at 500-mas resolution was done in the same configuration and 
season as Patnaik et al. (1992) or in the same configuration and contemporaneous with Browne 
et al. (1998). Those authors used phase calibrators 10° or less from their targets, cycled back to 
a calibrator every 6 minutes, and measured a 1-D positional accuracy of about 10 mas rms. The 
present work involves a calibrator 4° from the target (Figure 1) but observed with a 10-minute 
cycle time, so 20 mas was adoted as an estimate for the 1-D accuracy of the absolute VLA 
astrometry. In the VLBA case, the accuracy of the absolute astrometry was derived from analysis 
of the astrometric check source OQ208, as referenced to J1419+27. As Figure 1 shows, OQ208 
and J1419-I-27 are separated by about 3°, while the separation between J1419+27 and NGC 5548 
is about 2°. Thus the position error for OQ208 as referenced to J1419-I-27, plotted in Figure 5, 
serves as a cautious upper limit to the position error for NGC 5548 as referenced to J1419+27. 
Indeed, residual position errors of order 1 mas and 0.7 mas are expected for separations of 3° 
and 2°, respectively, due to known deficiencies in the 1980 lAU nutation model adopted at the 
correlator ( [Ma et al. 19981 ). 



3. Photometric Analysis 

Photometric monitoring results from Table 1 are plotted in Figure 6. Each figure panel shows 
results for a given angular resolution within 100 days centered on the mean observation epoch 
indicated. Data at 8.4 GHz are plotted as dark symbols, while data at 4.9 GHz appear as light 
symbols. The figure legend indicates which symbols encode which matched angular resolutions. 
Photometric comparisons should only be made among data at matched resolutions. To aid such 
comparisons the pairs of panels at the top, middle, and bottom of Figure 6 show data at matched 
resolution, with the left columns having the earlier mean epoch of observation. 

Figure 6 demonstrates that the radio nucleus of NGC 5548 is photometrically variable, on 
timescales that are both short (comparing within a panel) and long (comparing across paired 
panels). Specifically, the maximum flux densities are higher than minimum flux densities, by 



33 ± 5% between VLA observations separated by 41 days (panel [d]) and by 52 it 5% between 
VLA observations separated by 4.1 yr (panel [a] versus [b]). Panel (d) also establishes that the 
41-day "flare" is milder (19 it 5%) at 4.9 GHz than at 8.4 GHz. Furthermore, the spectral index 
a {S (X i^+") is a ~ -0.2 ± 0.1 prior to the 41-day "flare" but a ~ -hO.3 ± 0.1 during the "flare". 
Thus the "flare" is optically thick and milder at the lower frequency, traits typical for variable 
active nuclei. The "flare" is too poorly sampled to derive a meaningful variability timescale. The 
flux densities plotted in Figure 6 correspond to 8.4-GHz powers P^, ~ 2.5 — 8.3 x 10^^ W Hz~^, 
while the associated luminosities are Ljy = uP^^ ~ 2.1 — 7.0 x 10'^^ erg s~^. 

The radio sampling in Figure 6 is too coarse for useful comparisons with concurrent 
continuum monitoring from the ground or space. However, some relevant comparison can be made 
in connection with the 41-day radio variability between 1993 March 15 and April 25. First, that 
radio variability is much slower then the 1-day and 3-day variability reported for the hard and soft 
x-rays, respectively, during the period 1992 December to 1993 January ( Done et al. 1995| ). Second, 



the 33% change during the 41-day radio "flare" is comparable to the 40% changes documented for 
the continuum at 5100A during the period 1993 March 19 to May 27 ( |Korista et al. 1995| ). 



3.1. Implications for Black Hole Models 

Radio continuum emission can arise from (1) an advection-dominated accretion flow (ADAF) 
onto a black hole ( [Yi &: Boughn 1998| ), (2) jets emerging from an ADAF ( [Yi fc Boughn 1999| ), or 
(3) jets emerging from a standard accretion disk which itself emits no radio photons ( [Falcke fc| 
[Biermann 1996| ). Comparisions with these models will assume a black hole mass ~ 7 x lO"^ Mq 
( Peterson Sz Wandel 1999| ) and will try to match the spectral indices a, luminosities L^,, and 



photometric variability observed for the radio nucleus of NGC 5548, plus be consistent with VLBA 
emission being predominantly unresolved. 

The radio continuum from an ADAF is synchrotron emission fron an X-ray emitting and 
optically-thin plasma ( Yi fc Boughn 1998| , Yi &: Boughn 1999 ), and is characterized by three 



key traits. First, Equation 1 of Yi & Boughn (1999) predicts that the radio power should have 
an inverted spectral index a = -1-0.2, consistent with the index measured in the "flaring" state 
(a ~ -1-0.3 lb 0.1) but not in the the quiescent state {a ~ —0.2 it 0.1). Second, Equation 2.7 of Yi 
& Boughn (1998) implies that the 8.4-GHz photons should originate about 6 Schwarzschild radii 
from the black hole, placing them at radii less than 0.1 /xas, which is consistent with the VLBA 
emission being predominantly unresolved. Finally, Equations 1 and 3 of Yi & Boughn (1999) 
predict an 8.4-GHz luminosity Li, = uPi, ~ 7 x 10^^ erg s~^. This is about 30 — 100 times weaker 
than observed and argues for an alternative source of radio photons. Radio emission from jets is a 
plausible alternative, but before turning to that topic the existing evidence for outflow from the 
nucleus of NGC 5548 will be described. 

Figure 7 displays a VLA image at 8.4 GHz with resolution 660 mas = 350 pc and rms 



sensitivity ai = 12 /zJy ba^^. This deep image was obtained by adding the 1988-1989 VLA data 
lacking evidence for photometric variabihty (see panels [a] and [c] in Figure 6) and two further 
1989 observations at lower resolution. Areal integration over the figure yields a total flux density 
of 8.5 ± 0.3 mJy, corresponding to a total power at 8.4 GHz of 1.2 x lO'^^ W Hz~^. Figure 7 
shows that bipolar lobes straddle the radio nucleus and span 12 arcsec = 6.4 kpc. The lobes are 
very faint at 8.4 GHz, with peak intensities < 160 /iJy ba~^. The lobe morphologies confirm 
previous hints from images at lower resolution (Figure 3, Wilson &: Ulvestad 198^ Baum et al 



1993 , Kukula et al. 199q , Nagar et al. 1999| ) that the southern lobe is amorphous and the northern 



lobe is edge-brightened, possibly due to a bubble-like structure in 3-D. 

These synchrotron-emitting lobes could be driven by a bipolar wind emanating from the 



nucleus of NGC5548, as has been suggested for NGC3079 (Baan & Irwin 1995, Trotter et al. 



1998). Models for the absorption of the X-ray and ultraviolet continuum from NGC 5548 do 



imply outflowing ionized gas on subparsec scales ( [Mathur, Elvis, fc Wilkes 1995 ). The shape of 



the subparsec outflow (radial versus bipolar versus cylindrical) is not constrained, but its kinetic 
luminosity of 10^^ erg s~^ is four orders of magnitude higher than the photon luminosities of the 
radio lobes ( [Wilson fc Ulvestad 198^ ). If the subparsec outflow contains relativistic plasma, then 
Figure 4 could be used to trace the shape of the flow on parsec scales provided the reality of the 
apparent VLBA structure was verified (see Section 2). It is worth noting that weak synchrotron 
emission from relativistic plasma on scales 100-500 pc could be responsible for the source resolution 
effects, discussed in Section 2, that prompted the use of matched angular resolutions. 

Suppose the synchrotron luminosity of the jets on scales ^ 500 pc is powered by a rotating 
black hole accreting from a magnetized plasma. Then the jet luminosity can be estimated from 
the Blandford-Znajek process given values for the spin paramter a < 1 of the black hole, the 
efficiency ejet < 1 of the jet radio emission, and the mass accretion rate m, normalized to the 
Eddington rate ( |Yi &: Boughn 1999| ). Adopting m from an ADAF analysis of NGC 5548 leads to 
radio jet luminosities of Ly ~ 8.4 x lO^^a^ejet erg s~^. This process could easily match the required 
radio jet luminosites on scales ^ 500 pc, but the predictions are woefully imprecise compared to 
the predictions from the ADAF models. This process could also in principal match the kinetic 
luminosity of 10^^ erg s~^ for the subparsec outflow ( [Mathur, Elvis, fc Wilkes 1995| ). Alternatively, 
in the jet-disk symbiosis model of Falcke & Biermann (1996), the radio jet luminosity and the 
(standard) accretion disk luminosity are correlated. Gauging the latter for NGC 5548 by the 
nuclear bolometric luminosity of 5 x 10^^ erg s~^ ([Mathur, Elvis, &: Wilkes 1995 ), then the radio 



jet luminosity is as expected from that jet/disk correlation. Moreover, the spectral index of 
the quiescent radio emission is as expected for an inhomogeneous jet model, while the inverted 
spectrum during the 41-day "flare" could arise from a brightening in the optically-thick base of 
the jet. 



3.2. Implications for Star Cluster Models 

Terlevich et al. (1995) note the relevance of SN 1988Z to their models, as being an example 



of a (radio) supernova expanding into dense circumstellar material (Stathakis & Sadler 1991). It 



is thus of interest to compare the radio properties of this most distant known radio supernova 



(|Van Dyk et al. 1993D , with those of the radio nucleus in NGC 5548. Van Dyk et al. (1993) 
conclude that SN 1988Z is a rare, peculiar Type II supernova, which represents the endpoint of 
stellar evolution for a massive progenitor (~ 20 — 30 Mq) that underwent a period of high mass 
loss {'^ 10""^ Mq yr~^) before explosion. SN 1988Z achieved a maximum power at 8.4 GHz of 
about 2 X 10^^ W Hz~^, similar to the minimum power recorded in Figure 6 for the radio nucleus 
of NGC 5548. Van Dyk et al. (1993) present models quantifying the photometric evoluton of 
SN 1988Z at 8 and 5 GHz. (While the evolution appears to be smooth, the defining light curves 
are not well sampled in time.) Their models are plotted in Figure 8 with an arbitrary epoch shift 
for comparison with the VLA data from panel (d) of Figure 6. The radio supernova models appear 
able to account for the photometric variability of NGC 5548 observed on time scales of years. 
Those models cannot, however, account for the variability observed on time scales of tens of days. 



4. Astrometric Analysis 

All astrometric monitoring results from Table 1 are plotted in Figure 9. Each panel shows 
data at the matched angular resolution indicated, with the panels ordered as in Figure 6 to aid 
astrometric comparisons among data at matched resolutions. The origin for each panel is taken 
to be the position of the nucleus measured from the VLBA image at epoch 1995.337. Data at 
8.4 GHz are plotted as dark symbols, while data at 4.9 GHz appear as light symbols. The symbol 
shapes encode the matched angular resolution and have the same meaning as in Figure 6. 

Figure 9 implies that the radio nucleus of NGC 5548 is astrometrically stable, on timescales 
that are both short (comparing within a panel) and long (comparing across paired panels). Panels 
(a) and (b) demonstrate that the nucleus is astrometrically stable at 8.4 GHz, to an accuracy of 
20\/2 mas = 28 mas = 15 pc between VLA observations spanning both tens of days and 4.1 yr. 
Panels (c) and (d) verify this stability but with a poorer accuracy of 25\/2 mas = 35 mas = 19 pc. 
The nucleus is as astrometrically stable at 4.9 GHz as it is at 8.4 GHz, based on comparisons 
within panels (b) and (d). Finally, panels (e) and (f) imply that the nucleus is astrometrically 
stable at 8.4 GHz, to an accuracy of 1.3\/2 mas = 1.8 mas = 0.95 pc between VLBA observations 
spanning 13 days and 3.1 yr. 



4.1. Implications for Black Hole Models 

The observed astrometric stability is consistent with any models based on black holes, 
including ADAF and jet models, but an ADAF model was ruled out photometrically in Section 3.1. 
If the reality of the apparent VLBA structure in Figure 4 could be verified, then those images could 
be used to constrain the evolution of the lengths of, or features within, any jets or counterjets. 
Furthermore, if a VLBA image was acquired during a radio flare, then the peak in that image 
would arise from the optically-thick base of the jet and would be expected to be located closer to 
the black hole. 



4.2. Implications for Star Cluster Models 

Suppose that the blue luminosity from the Seyfert "nucleus" of the NGC 5548 arises from a 
compact star cluster. Then the absolute magnitude of the cluster can be equated to the historical 
minimum of Mb = —20.6 mag recorded for the Seyfert nucleus ( [Aretxaga &: Terlevich 1993| ). A 
star cluster with that magnitude will be characterized by an optical supernova rate of 0.6 yr~^ 
(Aretxaga & Terlevich 1994) and an effective radius r^, = 22 pc = 42 mas ([Young 1976, Melnick, 



Aretxaga, & Terlevich 1997). Each new optical supernova event should yield a radio counterpart 



if SN 1988Z is indeed as relevant to the starburst models as claimed by Terlevich et al. (1995). 

A circle with the effective radius r^ of the star cluster is shown in panels (a)-(d) in Figure 9, 
with the circle arbitrarily centered at the origin. Figure 9 implies that radio nucleus of NGC 5548 
does not move with the characteristic step size r^ between supernova events, as predicted by 
Melnick et al. (1997). The strongest constraints follow from the VLBA- VLBA astrometric stability 
at the 1-pc level for observations separated by 3.1 yr. However, in the star cluster model it is the 
massive stars which are the supernova progenitors, so if the initial mass function is biased toward 
massive stars and/or mass segregation is occuring, then the radio supernova would occur within a 
region of radius much smaller than rg = 22 pc = 42 mas ( [Melnick, Aretxaga, fc Terlevich 1997| ). 

The optical supernova rate of 0.6 yr"^ for NGC 5548 is moderately high, leading to a new 
supernova event about every 1.7 yr. This motivates a deep search for secondary radio sources at 
the matched VLBA resolution. Data from VLBA epochs 1998.445 and 1998.481 were combined 
to yield an image with rms aj = 41 /xJy ba~^, corresponding to about 100 times the strength of 



CasseopiaA were it placed at the distance of NGC 5548 (|Baars et al. 1977 ). No secondary sources 



stronger than 5c7/, or about 200 times CasseopiaA, were found in about 7600 VLBA beam areas 
within an annular region, centered as for Figure 9 and with inner radius 10 mas = 5.3 pc and 
outer radius 84 mas = 45 pc = 2re. This outer radius is well within the limitations to the field 
of view given the VLBA observing strategies ( Wrobel 1995[ ). The absence of secondary VLBA 



sources, in conjuction with the VLBA constraints on the astrometric stability of the radio source, 
requires the presence of no more than one radio supernova during a 3.1-yr interval. This condition 
is improbable if the detected VLBA source really is a radio supernova whose slow decline, based 



on the VLBA photometry, is shared by other radio supernovae. Thus either the radio supernova 
rate is lower than 0.6 yr~^ and/or, as mentioned above, the radio supernovae occur within a region 
of radius much smaller than rg = 22 pc = 42 mas. 



5. Conclusions and Future 

The VLA and VLBA were used to monitor the radio continuum counterpart to the optical 
broad line region (BLR) in the Seyfert galaxy NGC 5548. The nucleus is photometrically variable 
at 8.4 GHz by 33 ± 5% and 52 it 5% between VLA observations separated by 41 days and 4.1 yr, 
respectively. The 41-day changes are less pronounced (19 it 5%) at 4.9 GHz and exhibit an inverted 
spectrum (a ~ +0.3 it 0.1, S oc i'^") from 4.9 to 8.4 GHz. The nucleus is astrometrically stable 
at 8.4 GHz, to an accuracy of 15 pc between VLA observations separated by 4.1 yr and to an 
accuracy of 0.95 pc between VLBA observations separated by 3.1 yr. 

The photometric variability and astrometric stability of the radio nucleus is not consistent 
with star cluster models that treat the BLR as a compact supernova remnant and, for NGC 5548, 
required a new supernova event like SN1988Z every 1.7 yr within an effective radius rg = 22 pc. 
The radio supernova models for SN 1988Z cannot account for the VLA variability observed on 
time scales of tens of days. The combined VLBA evidence for astrometric stability, no secondary 
sources, and slow photometric evolution of the detected source mean that the radio supernova 
rate is lower than predicted and/or the radio supernovae occur within a region much smaller than 
predicted. Additional VLBA monitoring would thus be benefical, to provide further astrometric 
and photometric constraints on the detected radio source and to continue searches for secondary 
sources. 

The photometric variability and astrometric stability of the radio nucleus in NGC 5548 is 
consistent with models of jets on scales ^ 500 pc emerging from a black hole. If the subparsec 
outflow inferred from X-ray/UV absorber models contains relativistic plasma, then the VLBA 
images could be used to trace the shape of the flow on parsec scales. However, the reality of the 
apparent VLBA structure must first be verified. Tracing flow kinematics to 100 pc should also 
be feasible through an HST study of the narrow optical emission lines, because those lines are so 
dominated by gas within 70 pc = 130 mas from the nucleus ( [Kraemer et al. 1998 and references 
therein). 



A deep image at 8.4 GHz with a resolution 350 pc was obtained by adding data from 
quiescent VLA observations. This image shows faint bipolar lobes straddling the radio nucleus 
and spanning 6.4 kpc. The northern lobe is edge-brightened, suggesting a bubble-like structure in 
3-D. These lobes are plausibly driven by twin jets or a bipolar wind from the Seyfert 1 nucleus. 
The kinematics of the flow could be followed to kiloparsec scales if neutral or molecular gas was 
detected in absorption against the radio continuum lobes. Sensitive polarimetry of the radio lobes 
could also constrain the properties of an enveloping magnetoionic medium and the topology of the 



10 



magnetic field in a synchrotron-emitting bubble. 
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Fig. 1. — Phase reference sources near NGC 5548. 

Fig. 2. — Montage of VLA images of Stokes / emission from NGC 5548 at a matched resolution of 
500 mas = 260 pc (hatched circle). Contours are at intervals of 25 times the bottom contour level 
{BCL). Negative contours are dashed and positive ones are solid. Image peaks appear in Table 1. 
Upper row: Frequency 8.4 GHz. BCL = 100 ^Jy ba'^ Epochs left to right are 1988.913, 1988.989, 
1989.036, and 1989.077. Middle row: Frequency 8.4 GHz. BCL = 150 /xJy ba'^ Epochs left to 
right are 1992.814, 1992.858, 1992.891, 992.997, and 1993.025. Bottom row: Frequency 4.9 GHz. 
BCL = 150 /uJy ba~^. Epochs left to right are as for the middle row. 

Fig. 3. — Montage of VLA images of Stokes / emission from NGC 5548 at a matched resolution of 
1210 mas = 640 pc (hatched circle). Contours are at intervals of 22 times the bottom contour level 
{BCL). Negative contours are dashed and positive ones are solid. Image peaks appear in Table 1. 
Upper row: Frequency 8.4 GHz. BCL = 100 ^Jy ba~^ Epochs left to right are 1989.184, 1989.230, 
1989.279, and 1989.318. Middle row: Frequency 8.4 GHz. BCL = 150 ^Jy ba"^ Epochs left to 
right are 1993.162, 1993.203, 1993.249, 1993.277, and 1993.315. Bottom row: Frequency 4.9 GHz. 
BCL = 150 /uJy ba^^. Epochs left to right are as for the middle row. 

Fig. 4. — Montage of VLBA images of Stokes I emission from NGC 5548 at a frequency of 8.4 GHz 
and a matched resolution of 2.3 mas x 1.1 mas = 1.2 pc x 0.58 pc (hatched ellipse). Contours are 
at intervals of 22 times the bottom contour level of 150 /^Jy ba~^. Negative contours are dashed 
and positive ones are solid. Total flux densities appear in Table 1. Left: Peak 1.93 mJy ba~^. 
Epoch 1995.337. Middle: Peak 1.09 mJy ba'^ Epoch 1998.445. Right: Peak 0.73 mJy ba^^ 
Epoch 1998.481. 

Fig. 5.— Observed minus expected VLBA position for OQ 208 at UT epochs 1995.337, 1998.445, 
and 1998.481. 

Fig. 6. — Photometric monitoring of the radio continuum nucleus of NGC 5548 at the matched 
angular resolutions and mean epochs indicated. Data are from Table 1. Small circles and squares 
symbolize peak flux densities at a resolution of 500 mas = 260 pc (see Figure 2). Big circles and 
squares indicate peak flux densities at a resolution of 1210 mas = 640 pc (see Figure 3). Small 
diamonds correspond to integrated flux densities at a resolution of 2.3 mas = 1.2 pc times 1.1 mas 
= 0.58 pc (see Figure 4). 

Fig. 7. — Deep VLA image of Stokes / emission from NGC 5548 at a frequency of 8.4 GHz, resolution 
of 660 mas = 350 pc (hatched circle), and a mean epoch of 1989.192. Contour levels are at intervals 
of 22 times the bottom contour level, which is 30 jiiy ba~^. Negative contours are dashed and 
positive ones are solid. Grey scale ranges linearly from 30 to 180 fxiy ba~^. 

Fig. 8. — Photometric monitoring of the radio continuum nucleus of NGC 5548 during the 41-day 
"flare" , for comparison with models for the photometric evolution of the radio supernova SN 1988Z. 
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Fig. 9. — Astrometric monitoring of the radio continuum nucleus of NGC 5548 at the matched 
angular resolutions and mean epochs indicated. Data are from Table 1. Symbols have the same 
meaning as in Figure 6. Section 4.2 describes the significance for the star cluster models of the 
circle in panels (a)-(d). 
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Table 1. Photometric and Astrometric Monitoring 
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Epoch 
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Freq. 
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(jiiy ba 
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q;j2ooo 
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8.4 


.279 .. 


. VLA 


1210 
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4.9 


.891 .. 


VLA 


500 


8.4 


.891 .. 


. VLA 


500 
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"Peak in mjy ba ^. 

''Total in mJy from areal integration. 



